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The modeling and analysis of the dynamics of complex systems often requires to employ non-Markovian 
stochastic processes. While there is a clear and well-established mathematical definition for 
non-Markovianity in the case of classical systems, the extension to the quantum regime recently caused a 
vivid debate, leading to many different proposals for the characterization and quantification of memory 
effects in the dynamics of open quantum systems. Here, we derive a mathematical representation for the 
non-Markovianity measure based on the exchange of information between the open system and its 
environment, which reveals the locality and universality of non-Markovianity in the quantum state space 
and substantially simplifies its numerical and experimental determination. We further illustrate the 
application of this representation by means of an all- optical experiment which allows the measurement of 
the degree of memory effects in a photonic quantum process with high accuracy. 

In recent years the problem of characterizing non-Markovian dynamics in the quantum regime has initiated an 
intense debate. A series of diverse definitions along with measures of quantum memory effects have been 
proposed, invoking many different mathematical and physical concepts and techniques. Examples are char- 
acterizations of non-Markovianity in terms of deviations from a Lindblad semigroup\ of the divisibility of the 
dynamical map^, of the dynamics of entanglement^ and correlations^ with an ancilla system, and of the Fisher 
information^. 

In this work we focus on the measure of non-Markovianity introduced in Refs. [5, 6] which characterizes non- 
Markovianity in terms of the information exchanged between an open system and its environment. To quantif)^ 

this information one uses the trace distance V{p^,p2) = - Tr — P2 1 between two quantum states represented by 

density matrices pi and p2, which provides a measure for the distinguishability of these states^"^. The dynamics of 
an open quantum system can be described formally by means of a dynamical map which maps any initial state 
p(0) to its time-evolved state p(t) = ^tp(O) at time f^. The time evolution over some (finite or infinite) time 
interval [0, T] is then described by the one-parameter family O = {O^ | 0 < t < T, = 1} of dynamical maps. A 
quantum process given by such a family of dynamical maps is said to be Markovian if for all pairs of initial states 
the trace distance V[p^(t) ,P2 (0 ) ^ monotonically decreasing function of time. In physical terms this means that 
the distinguishability between the states decreases monotonically which can be interpreted as a continuous flow of 
information from the open system to the environment. Conversely, a quantum process is non-Markovian if there 
is an initial pair of states for which the trace distance increases over certain time intervals. During these time 
intervals the distinguishability of the states thus increases, which corresponds to a flow of information from the 
environment back into the open system signifying quantum memory effects. 

The above characterization of quantum non-Markovianity leads to a measure for the degree of non- 
Markovianity which is defined by 



A/'(0) = max 



/ dt G{t,p^,P2), 

J (7>0 



(1) 



where 



0{t,p,,P2)=jV{^t{Pl)MPl)) 



(2) 



SCIENTIFIC REPORT; | 4 : 6327 | DOI: 1 0.1 038/srep06327 



1 



denotes the time derivative of the trace distance between the pair of 
states at time t. In Eq. (1) the time integral is extended over all time 
intervals in which this derivative is positive, and the maximum is 
taken over all pairs of orthogonal initial states pi ± p2- Thus, the 
measure accumulates the total increase of the distinguishability 
occurring in the time evolution of two initially orthogonal quantum 
states. We recall that two quantum states pi and p2 are said to be 
orthogonal if their supports, i.e. the subspaces spanned by their 
nonzero eigenvalues, are orthogonal which is equivalent to 
^(Pi^Pi) = 1- This impUes that optimal state pairs exhibiting a max- 
imal backflow of information during their time evolution are initially 
distinguishable with certainty, and thus represent a maximal initial 
information content. 

This measure for non-Markovianity was originally defined in Ref 
[5] in terms of a maximization over all pairs of quantum states. 
However, as demonstrated in Ref. [11] the maximization can be 
restricted to pairs of orthogonal initial states. Although this result 
greatly simplifies the mathematical representation of the non- 
Markovianity measure, its determination still requires the maximiza- 
tion over pairs of quantum states. Here, we derive a much simpler 
representation for the measure which is particularly relevant for its 
experimental realization since it only requires a local maximization 
over single quantum states, the second state being an arbitrary fixed 
reference state taken from the interior of the state space. This rep- 
resentation will further be employed in an all -optical experimental 
setup for the measurement of the non-Markovianity of a photonic 
quantum process. 

Results 

Local representation of quantum non-Markovianity. In the 

following we denote the Hilbert space of the open quantum 
system by H (assumed to be finite -dimensional to simplify the 
presentation) and the corresponding state space, given by the 
convex set of density matrices with unit trace, by S{H). To 

formulate our main theoretical result we define S{H) to be the 
interior of the state space, i.e. the set of all quantum states po for 
which there is an g > 0 such that all Hermitian operators p with 
unit trace satisfying V{p,Pq)<s belong to S{H). We further 
define 



So{n) = {A\A^O,A = A\ Tr A = 0} 



(3) 



to be the set of all nonzero, Hermitian and traceless operators on 

7Y. Considering any fixed reference state Po eS {Ti) we can now 
introduce a particular class of subsets of the state space: A set 
dU^Poj^SiH) not containing po is called an enclosing surface 
of Po if ^nd only if for any operator A e So{7{) there exists a real 
number A > 0 such that 



Po + ^A G dU{pQ 



(4) 



Note that by definition po itself is not contained in dU{po) and 
that the full set dU(po) is part of the state space. It can be easily 
seen that any state from the interior of the state space has an 
enclosing surface. For example, since po is an interior point of 
the state space there is an e > 0 such that the set of states p 
defined by V{p,Pq)=s represents a spherical enclosing surface 
with center po. However, an enclosing surface dU(po) can have 
an arbitrary geometrical shape, the only requirement being that it 
encloses the reference state in all directions of state space. An 
example is shown in Fig. 1(a). Using these definitions, we can 
now state our central result. 

Theorem. Let pg e 'S{H) be any fixed state of the interior of the state 
space and dU(po) an arbitrary enclosing surface of po- For any 
dynamical process ^, the measure for quantum non-Markovianity 
defined by Eq. (1) is then given by 




(a) 



(b) 



Figure 1 | Illustration of an enclosing surface dU{pQ) (a) and of a 
hemispherical enclosing surface dU{pQ) with disconnected boundary (b) 
for an interior point pQ of the state space S{H). 



where 



J\f{^)= max / (it (7(t,p,Po), 
|P(cD,(p),0,(Po)) 



a{t,p,pQ)=' 



^{P.Po) 



(5) 



(6) 



is the derivative of the trace distance at time t divided by the initial 
trace distance. 

Proof. Let p e dU(po). Applying the Jordan-Hahn decomposition^ 
to the operator p — Po one concludes that there exists an orthogonal 
pair of states pi and p2 such that 

P-Po 



P\-P2 = 



and, hence, we have 



D(<I>,(p,)A(P2)) = 



(7) 



(8) 



by the linearity of the dynamical maps and the homogeneity of the 
trace distance. This shows that G(t,Pi,p2) = a-(t,p,Po). It follows that 
the right-hand side of Eq. (5) is smaller than or equal to A/'(0) as 
defined by Eq. (1). Conversely, suppose Pi, Pi are two orthogonal 
states. Since Pi— Pi ^ ^o(^)j there exists > 0 such that p = Po + 
KP\ ~ Pi) ^ ^U(po), by definition of an enclosing surface. Thus, 
one obtains Pi — P2 = (p ~ Po)/^- Since pi ± P2 we find 
V{p,Pq)/X = V{pi,P2) = 1 and, hence, X = V[p,Pq). Thus, we are 
again led to Eq. (7) and to a{t,Pi,p2) = o-(t,p,Po). This shows that 
the measure A/'(^) as defined by Eq. ( 1 ) is smaller than or equal to the 
right-hand side of Eq. (5) which thus concludes the proof 

The statement of the theorem can be easily understood for the case 
of a qubit, representing the corresponding state space by means of the 
Bloch ball. The representation (1) requires to perform a maximiza- 
tion over all pairs of orthogonal quantum states, i.e. over all pairs of 
antipodal points on the surface of the Bloch ball. The theorem states 
that this maximization can also be carried out in the following equi- 
valent way. Choose a fixed point po in the interior of the Bloch ball 
and a small surface ^L^(po) which encloses po, and take the maximum 
over all points p G dU{po). To see the equivalence, consider the 
straight line joining po with p. This line can be moved by parallel 
translation (without changing the trace distance) in such a way that it 
is centered at the origin of the Bloch ball (maximally mixed state). 
Stretching the line by an appropriate scaling factor (which is equal to 
the inverse of the trace distance between po and p) its endpoints 
become antipodal points of the Bloch ball. Conversely, any pair of 
antipodal points on the surface of the Bloch ball corresponds in this 
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HWP QWP FP IF attenuator QP BS PBS mirror SPD 

Figure 2 | Experimental setup. Key to the components: HWP - half-wave 
plate, QWP - quarter- wave plate, FP - Fabry-Perot cavity, IF - interference 
filter, QP - quartz plate, (P)BS - (polarizing) beamsplitter, SPD - single 
photon detector. 

way to a point p belonging to dU{po), which leads to the theorem for 
the particular case of a qubit. 

Experimental realization. We have used Eq. (5) to develop an 
experimental scheme for the determination of the degree of non- 
Markovianity in a photonic quantum process. In the experiment the 
open quantum system is given by the polarization degree of freedom of 
a single photon coupled to the frequency degree of freedom 
representing the environment. While in many system-environment 
models the system and its environment are represented by different 
physical entities, the present situation in which system and 
environment are formed by different degrees of freedom of the same 
particle is not uncommon in the theory of open system. A typical 
example in this context are experiments on trapped ions where also 
the internal electronic degree of freedom couples to the motional 
degree of freedom of the ion^^. In our experiment the decoherence of 
the polarization degree of freedom is due to birefringent quartz plates 
in the optical path of the photon which induce a coupling between the 
polarization and frequency degrees of freedom and lead to dephasing 
of superpositions of vertical and horizontal polarization states. This 
dephasing strongly depends on the structure of the frequency 
spectrum which can be efficiently controlled by the tilt angle of a 
Fabry-Perot cavity, producing a bimodal spectrum. Hence, the non- 




(a) 



(b) 



(c) 




Figure 3 | Experimental results for the increase of the trace distance 
between 17 5 A and 318y2 for = 0.64 for states on the enclosing surface of 
reference state pj (a), pi (b) and pairs of orthogonal states (c). The 

corresponding 0ioc- averaged increase with respect to local spherical 
coordinates is shown in (d), (e) and (f). Error bars show the standard 
deviations. 



Figure 4 | The same as Fig. 3 for = 0.22. 

Markovianity in our experiment is due to the presence of a structured 
environment which is a typical cause for memory effects. 

The experimental setup is depicted in Fig. 2. With the help of a 
frequency doubler a mode-locked Ti:sapphire laser (central wavelength 
780 nm) is used to pump two 1 mm thick BBO crystals to gene- 
rate the maximally entangled two-qubit state {\H ,V) — \V ,H)) / Vl 
with \H) and |V) denoting the horizontal and vertical polarization 
states, respectively^^. A fused silica plate (0.1 mm thick and coated with 
a partial reflecting coating, with approximately 80% reflectivity at 
780 nm) serves as a Fabry-Perot cavity (FP) which in addition can 
be tilted to generate different dynamical behavior^^. The cavity and a 
consecutively placed interference filter (IF) (FWHM about 3 nm) sin- 
gle out two peaks near 780 nm of width g = 7.7 X 10^^ Hz each which 
are separated by Aco = 7.2 X 10^^ Hz. The relative amplitude A^, of the 
two peaks depends strongly on the tilt angle a whereas the other 
quantities are almost constant. A polarizing beamsplitter (PBS) 
together with a half-wave plate (HWP) and a quarter-wave plate 
(QWP) are used as a photon state analyzer^^. 

Photon 1 is directly detected in a single photon detector at the end of 
arm 1 as a trigger for photon 2. The optical setup in part a, b and c (see 
Fig. 2) is used to prepare arbitrary quantum states of photon 2 needed 
for the sampling process^^. This set-up conveniently allows to prepare 
any single pure photon polarization state (in arm 2c) and reference 
states {2a along with 2b) together with arbitrary enclosing surfaces 
which can be controlled by changing the relative amplitudes of the 
attenuators built in in each arm. The path difference between each arm 
is about 25 mm to ensure that the mixture of the three parts is classical. 

After the preparation photon 2 passes through birefringent quartz 
plates of variable thickness which couple the polarization and fre- 
quency degree of freedom and lead to the decoherence of superposi- 
tions of polarization states. The birefringence is given by An = 8.9 X 
10"^ at 780 nm. The thickness of the quartz plates simulating differ- 
ent evolution times ranges from 75k to 318>^ in units of the central 
wavelength of the FP cavity. 

Employing the Bloch vector representation, the set of polarization 
states can be conveniently parametrized by means of spherical coor- 



nTTirrn'' | 




(a) 



(b) 



(c) 
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Figure 5 | The same as Fig. 3 for = 0.01. 
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Table 1 The quantum 


non-Markovianity 


measure for the three dynamics 


obtained from the experimental data 


in comparison to the 


theoretical value 










A. 


ATtheo 








0.64 


0.59 


0.59 ± 0.01 


0.59 ± 0.02 


0.59 ± 0.02 


0.22 


0.21 


0.21 ± 0.01 


0.21 ± 0.02 


0.21 ± 0.02 


0.01 


0 


0.001 ± 0.013 


-0.005 ± 0.008 


0.0002 ± 0.0015 



dinates r = (r, 0, (j)). We apply the local representation to two ref- 
erence states to determine experimentally the degree of non- 
Markovianity for three dynamics characterized by the relative ampli- 
tudes = 0.64, 0.22 and 0.01, ranging from non-Markovian to 
Markovian evolutions, and compare the results with the outcome 
for pairs of orthogonal initial states. The reference states pj and pi 
used in the experiment are given by 

ri= f 0.20, -71, —71 I, f 0.88, — 7i, — 71 V (9) 

° V 2 50 y ° V 50 50 y ^ ^ 

Reference state pj is thus located inside the equatorial plane, whereas 
the second reference state lies in the northern hemisphere close to the 
boundary. The enclosing surfaces are determined by the convex 
combination O.S-pJ' ^ + 0.7-p of the reference states and any pure 
state prepared in arm 2c. These sets thus contain only mixed states. 
We measured a total of 5000 states on the surface for each reference 
state which are characterized by the azimuthal and polar angles of the 
pure states. The associated angles 0 and (j) are located on a lattice with 
equal spacing of 27i/100. 

The outcomes of the measurements are presented in Figs. 3, 4 and 
5. The increase of the trace distance between 17 51 and 318/1 for any 
state on the enclosing surface for the two reference states is shown in 
Figs. 3(a)-5(a) and 3(b)-5(b) using color coding. Note, that the 
colored surfaces in these figures are non- spherical and not centered 
at the origin. By contrast, the ordinary Bloch spheres depicted in 
Figs. 3(c)-5(c) show the measurement outcomes for pairs of ortho- 
gonal initial states. 

Defining spherical coordinates rioc = (Hocj ^Iocj hoc) with respect 
to local coordinate systems centered at the position of the two ref- 
erence states pj and pg, one recovers the polar symmetry present for 
pairs of orthogonal states, see Figs. 3 (a) -(c). One may therefore 
average the outcomes over the polar angle (/>ioc along lines of latitude. 
To this end, we introduced an appropriate binning on the z-axis and 
determined the average increase which we then assigned to the azi- 
muthal angle ^loc associated to the mean z- value in the bin. In addi- 
tion, we allocated the standard deviation to each of the averaged 
outcomes. The resulting data are depicted in Figs. 3(d)-5(d) and 
3(e) -5(e) and show the same characteristics as the (^-averaged 
increase of pairs of orthogonal states displayed in Figs. 3(f) -5(f). 
Note that the directional dependence of the trace distance originating 
from its property of depending only on the difference of two states 
can be nicely seen for example in Figs. 3(d) -(f). 

The maximal increase of the trace distance for the two reference 
states obtained from the (^loc-averaged data as well as for pairs of 
orthogonal states are given in Tab. I. The experimentally determined 
values are in very good agreement with the predictions of the theor- 
etical modeP^, demonstrating the experimental feasibility and the 
accuracy of the method based on the local representation (5). 

Discussion 

The representation of quantum non-Markovianity given by Eq. (5) 
bears several important mathematical and physical consequences. 
First, it demonstrates that the non-Markovianity measure can be 
determined by maximization over single quantum states p taken 
from an arbitrary neighborhood of a fixed state po in the interior 



of the state space. Thus, Eq. (5) provides a local representation of 
non-Markovianity, showing that quantum memory effects can be 
detected locally by sampling single states from an arbitrary enclosing 
surface of a fixed reference state. Note that the theorem cannot be 

applied to infinite dimensional Hilbert spaces since S{7{) is empty in 
this case. 

Second, the choice of the fixed reference state po is completely 
arbitrary, the only condition being that it belongs to the interior of 
the state space. Thus, the non-Markovianity of a dynamical process is 
indeed a universal feature which appears everywhere in state space: 
The information about non-Markovian behavior is contained in any 
part of the state space which supports the intuitive idea that quantum 
memory effects represent an intrinsic property of the dynamical pro- 
cess. This fact is particularly relevant when dealing with a dynamical 
process that has an invariant state in the interior of the state space. It 
is then of great advantage to choose po as this invariant state such that 
only the sampled states p e dU{po) evolve nontrivially in time. 

Third, the theorem greatly simplifies the analytical, numerical or 
experimental determination of the non-Markovianity measure. In 
particular, it shows that it is not necessary to scan the whole state 
space in order to find an optimal pair of quantum states, but rather to 
sample the states of an enclosing surface of a fixed interior point of the 
state space. From the proof of the theorem we also see that it suffices if 
the enclosing surface contains all directions emanating from the fixed 
reference state po exactly once, i.e. if Eq. (4) holds for exactly one >^ > 0. 
It is even sufficient if this equation holds for either A or —A. Therefore, 
the theorem is also valid if dU{po) is replaced by a hemispherical 
enclosing surface dll^Po) which we defme as follows. A set 
dll^Po) (^SiTL) is said to be a hemispherical enclosing surface of po 
if and only if for any Ae 8q{1-L) there exists exactly one real number 
> 0 such that either p^ + kAe dU^Po) or p^ — kAe dU^Po). A hemi- 
spherical enclosing surface thus contains all directions, given by opera- 
tors A E So{?{), only once. Moreover, it needs neither be smooth nor 
connected (see Fig. 1(b) for an example) which makes this character- 
ization particularly useful for noisy experiments. 

Summarizing, we have developed a representation of the measure 
for quantum non-Markovianity which fully reveals the locality and 
universality of memory effects in the quantum state space. These prop- 
erties are clearly reflected in the results of our photonic experiment. 
The experiment illustrates that the measure can be obtained efficiently 
in an arbitrary neighborhood of any fixed state in the interior of the 
state space, that its determination only requires a maximization over a 
single input state, and that optimal quantum states featuring maximal 
backflow of information can always be represented by mixed states. 
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